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ON TRIPLY ORTHOGONAL CONGRUENCES* 

BT 

JAMES BYRNIE SHAW 

Introduction 

1. A vector field is defined by the value of a vector a at each point of the 
field. This vector a is a function of p, the vector to the point from some 
given origin. If <r is taken as tangent to a curve at each of the points of the 
space considered, these tangents will envelop a congruence of curves, the 
vector lines of the field of a . The tensor or length of a is not determinable 
from the congruence, the congruence depending only upon the unit vector Ua . 
The congruence itself however determines many properties of the field, which 
are of course purely geometric properties, and in physical phenomena are not 
dependent on physical facts except through the law by which Ua is determined. 
It is purposed to consider some of these properties which belong to the con- 
gruences determined by three unit vectors which are everywhere mutually 
perpendicular. These will generally be indicated by a, j8, y; but in case 
they are the tangent, principal normal, and binormal of a curve, we will use v 
with a proper subscript to indicate the principal normal, and n with the same 
subscript to indicate the binormal. For instance, if we are discussing a curve 
of the a congruence the normal will be v^ and the binormal /*„ • 

The algebra used throughout is Quaternions, all the usual symbols of that 
algebra being introduced with no explanation. However it should be noted 
that such an expression as afi is always a quaternion product. 

The symbol <l> and the symbol 6 are throughout linear vector operators, — 
in the case of related to certain vectors which are subscripts of ^. For 
instance we define for any vector ri :<}>,, = — S {) V • V • The occurrence of 
such operators is due to the fact that if the point under discussion at the end 
of p is displaced to p + dp, the vector rj becomes ri + drj where 

drj = — Sdp^ ■ ri = <l>^dp. 

In every instance V operates upon every variable that follows it unless it 
carries a subscript, in which case it operates upon everything with that sub- 
script whether preceding or following, the subscripts being dropped after the 
differentiation. 



* Presented to the Society, December 28, 1918. 
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2. We recall the Serret-Frenet formulas in this notation, indicating the 
curvature by c with proper subscript and the torsion by t with subscript. 

For instance, we have for a curve of the a congruence 

"a Oa. = 't>aOC, M« <« " ««„ = <f>y a , " "« <a = ^M « ' 

since <t>a « is the vector rate of change of « , dv/ds , in the direction of a , ^„ a 
the rate of change, dv/ds, of v, for displacements in the direction of a, etc. 
If we set a)„ = c„ /*„ + i„ a then a)„ is the rectifying line of the curve and we 
have (omitting the subscripts for brevity) 

dv dfj, da 

ds = ^'"'' d^=^'*"'' ^ = ^"«- 

In case we have rj = la + m^ + ny where I, m, n are constant and w retains 
its meaning. 

The system a, ;8, y 

3. Since a , ;8 , 7 are unit vectors Sada = = S^d^ = Sydy , or Sa<^„ dp = d 
for all dp, and similarly for the others. Hence we have for any vector, 
Sa<j>^ () =0,or^lo; = where ^' is the conjugate, or transverse operator, 
and likewise 4>'^^ = Q, ^'yy = 0. 

Since -Sa/3 = etc. we have Sad^ + S^da = , or 8a<i>^ ( ) + S/3^„ ( ) = , 
ov <j>'^a = — <j>'a^ with similar equations for the others. 

Since a^ = y we have V7 = Vi ai /3 + Vi 01^1 , and since for any two 
vectors X , /i we have always \p. = — /A + 2S\p, we obtain from the equation 
just written 

V7 = Vi ai i8 - Vi /3i a + 2Vi SajSi, 

or, remembering Vi Sa^i = — Vi Sai (8 , 

V7 = Vi«i/3- Vi/3ia + 2^;;8. 

We obtain the similar equations for the other vectors a , (8 by interchanging 
the symbols cyclically. This equation shows that the quaternion V7 is fully 
determined from a and j8. Taking the scalar and writing the expressions 
so that the subscripts may be omitted, the convergence of 7 is 

*SV7 = SfiVoi - 8aV^, or Sy [ VaVVa + V^Vv^ + F7FV7] . 

The expression in the bracket being the same for all three vectors, we will set 

2e = VaVVcc + V^Vv^ + F7FV7 
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and have the three convergences in the form 

<SVa = 2Sae, Sv^ = 2S^e, 2SVy = 28y€. 
Hence 

2e = — aSVa — /3SV/3 — ySVy • 
Since 

FaFVa = SaV • a — Vi Sai a = SaV • a — ^^Saa = — ^<, a = — c„ >»<, 

it becomes evident that e is half the sum of the three vector curvatures, with 
sign reversed, of the three curves for a, j8, 7 at the point considered. 

4. We have, by taking the vector of V7, the curl of 7: 

Vvy = FViaiiS - FVi 1810; + 2^1/3 = <^„;8 - <l>^ a + 2V€y . 

The last member is found by expanding the first two vectors by the usual qua- 
ternion formulas, introducing ^ , and substituting the values of Sv<x , Sv^ • 
This form will be further reduced. 

5. Multiplying V7 into 7 and taking the scalar we have (since S'K<I>' n 

= SfJKJA) 

SyVy = SaVa + S^V^ + 2S;8^„ 7 . 

If now we set 2p = SaVa + Sj8Vj8 + S7V7 we have 

Saij>fi 7 = — iSj8<^a T = P ~ S7V7 , S^<l>y a = — Sy<l>^ a = p — SaVa , 

S7^.j8 = - Sa<t>yp = p - Sj8V/3. 

Multiplying a into V7; taking the scalar, reducing by previous formulas, 
and writing the corresponding formulas cyclically, 

1S7 Va = SP<j>^ a = — Sa^^ a = c„ S^p^ , 

SaVP = Sy<t>pP = - Sp(j>y^ = c^ Syvp, 

S/3V7 = Saft>y 7 = — iS7^a T = c,, SaUy . 

Multiplying by j8 and taking the scalar, and writing corresponding formulas : 

SyV^ = S^<t>a^ = - Sa<t>p ^ = - c^ Sav^ , 

SaVy = Sy<t>^ 7 = — S^<t>y y = — Cy Sfipy , 

SfiV<x = Sa(t>y a = — 8y<t>a a = — c„ iS7»'„ . 

6. From these results it is easy to see that if p is the principal normal of 
the a curves and fi the binormal, c the curvature and t tlje torsion, p./^ referring 
to a, V , n, 

c = - SmV" = TVaVVoi, t = p,^ - SaVa, c" = - VaVVa, 

p = — UVaVVoi, n = — UVaVaVVoi, a)„ = p,^ a + Vv<x, 

(SaV)FVa 



Cfi = — VaVaVVce = aSaV<x + FVa , t = — SaVa — 8 



VaVVoc 
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Since the determination of the normal and the curvature depends only 
upon an expression of the f orm ^, rj for any unit vector rj , similar expressions 
for the curves of j8 , 7 , f , and /x may be written down at once except for t . 
The determination of t depends upon the differentiation of the unit normal in a 
direction perpendicular to it, and consequently cannot be found as simply 
in the general case. However in certain special cases to be studied, the 
determination can be simplified. This difficulty of determination is also 
evident when we notice that the expression for t contains p , which depends 
upon the curl of /3 and 7 as well as on that of a . In general the quantity p 
depends upon what ^ and 7 are, and in the formula for t it is understood that 
in finding p, ^, and 7 must be the normal and binormal respectively. 

In case ^ and 7 are not the normal and the binormal, let the angle from 
the normal to ;8 be w , then 

j8 = cos w • V -\- sin w • n, 
V/3 = cos w • Vf + sin 1^ • Vm + Vi wi • 7 , 
7 = —smw-v + cos w • n, 
V7 = — sin w • V + cos w • Vm ~ Vi wi • |8 . 
Hence we have 

«Si8V/3 4- S7V7 = SpVv + S/xVm - 2SaVw, 

and referred to a , /3 , 7 we have for the value of p in terms of p,^ (the value 
of p in terms oi a,v , n) 

2p = iSaVa + S^VP + /S7V7 = 2p^^ - 2SaVw, 

or p = py^ + derivative of w for displacements in direction a . That is to 
say the quantity p belonging to arbitrary axes will be the corresponding 
quantity for the intrinsic axes of the curve plus the angular rate of rotation 
of the arbitrary system (8, 7 about a, for displacements along a. If this 
angular rotation is zero, so that (3 maintains a fixed angle with the normal, 
then there is no change in p . The torsion of a in this case can be found from 
the a , /3 , 7 trihedral. 

In particular if we consider the (8 curves and let o; be the principal normal 
for the j8 curve and also for the 7 curve then for any ^ curve the torsion is 
p — Sj8V/3 , and likewise for any 7 curve the torsion is p — SyVy • It follows 
that the sum 

tp +ty = 2p - S^V^ - S7V7 = SaVoc 

and this is independent of the particular curves ^ and 7 provided a is their 
principal normal.* Hence the sum of the torsions of any two orthogonal 

* R. A. P. Rogers, Some differential 'properties of the orthogonal trajectories of a congruence 
of furves, mth an appUcaiion to airl and divergence of vectors, Proceedings of the 
Roj-al Irish Academy, 29A (1912), 92-117. 
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curves for which a is the principal normal is also Sa^a. For, the tangent 
to any curve of which a is the principal normal can be written j8 cos u + ysinu. 
Hence the normal a times the curvature c is given by 

— )S (j8 cos M + 7 sin m) V • (|3 cos m + 7 sin u) 

= a{cp cos^ u + Cy sitf u + sinu cos u) ( Sy<j>^ jS + S/3^„ 7) 

+ ;8 (sin M cos M (<S|8V« + SjS^y i3) + sin^ iiSyVu) 

— 7 (sin « cos u{SyVu — Sy^^ 7) + cos^ uS^Vu) . 

The terms in ^ and 7 must however vanish by hypothesis so that in as much as 
u is an arbitrary angle we must have along the curve in question Vm = 0, 
since 8^<t>y |3 = — Sy^p /3 = Sya = , and Sy<l>p y — — S^<t>y 7 = 0. Hence 
along any curve for which the vector a is the principal normal the tangent 
will make fixed angles with ^ and 7 , or what is the same thing the trihedral 
maintains a constant relation to the fundamental trihedral. Hence for any 
two such curves with perpendicular tangents, the sum of the torsions is SaV<x .* 

It follows that the quantity p is the sum of the mean torsions of curves 
normal to each of the three vectors a, jS, 7, respectively. These curves 
cannot be taken so as to serve for different principal normals, unless they are 
straight lines. For instance the two used for a as principal normal cannot be 
used also for /3 as principal normal. 

Further the curvature of a curve with a as principal normal is 

C(3 cos^ u + Cy sitf u + sin u cos « [ SjS ( 0„ + ^1) 7 ] . 

If we add to this the curvature of a perpendicular curve of the set we have 
Cp + Cy, which does not contain u, and is the same therefore for all such 
perpendicular pairs; it reduces to iS/3<^„j3 + Sy^^y, that is to — mi (</)„) 
where mi(<^) is the first scalar invariant of the operator </>, and for <t>a is 
+ 8'Voi ■ Hence the sum of the curvatures of the perpendicular pair of 
normal curves is — SVa .* 

Since SjSVa. is the projection of the vector curvature of a on 7, the per- 
pendicular to both a and /?, and Sy^a. the projection of this vector curvature 
on — j3, the common perpendicular to 7 and a, we see that the curl of a 
consists of a vector along a and a vector in the plane perpendicular to a , the 
latter being the numerical curvature times the unit vector which is the bi- 
normal. Hence the components of (8 along the unit normal and the unit 
binormal are — SyVa/c and — S^'^a/c, and the corresponding components 
of 7 are S^Va/c and — SyVcc/c. In case FaFVa = these components 

* R. A. P. Rogers, Some differential properties of the orthogonal trajectories of a congruence 
of curves, with an application to curl and divergence of vectors. Proceedings of the 
Royal Irish Academy, 29A (1912), 92-117. 
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become indeterminate. In case SaVa = , a is perpendicular to FVa and e 
becomes merely TV^a . Also we notice that in any case 

p = t^ + iSaVa + dw/ds . 

The operators <t>a, <i>^, 4>y 

7. Since ^^a = 0, the transverse of <j>a has at least one zero root, and 
therefore 0„ has at least one zero root. The form of <i>a must then be either 

(a) - mz + ysu, U 7^ 0, f2 5^ 0, VU f2 ?^ 0, 

Q>) |8' iSf, f ?^ 0, where jS' is perpendicular to a, j8'^ = — 1 , 

(c) 0. 

In case (c) a is constant and its congruence consists of parallel straight 
lines. The ;8 and y curves lie in parallel planes perpendicular to a, that is a 
is their common binormal. This is evident also from the Serret-Frenet 
formulas. For we have 

<^a=-S()V«=-'S()V^7 = V^<i>, - Vy<t>p , 

and if this vanishes its transverse vanishes, that is <j>'p Vy — ^^ F/3 = , whence 

^3 a = , <f>'ya = 0, and as 0^ j8 = , ^^7 = 0, 

<t>p = T-SI, <t>y= - |8S?, 4>pp = yS^^, <(>yy = - ^7, 

so that the binormals are a in each case. The j8 curves in their plane are 
respectively orthogonal to the 7 curves in the same plane. The curvature 
of the ;8 curves is the projection of J on their normal, and likewise for the 
curvature of the y curves. The projection of f on a is the quantity — p . 
In case (b) we find that for some vector t 

<i>p = — aS^ cos u + ySr , <t>y = — aS^ sin u — fiSr , 

where j8' = )3 cos m + 7 sin « . In this case da = + ^' S^dp for any dp . 
But if dp is taken along the a curve da = c„ v^ hence + j8' = »». , a/3 = |x„ , 
and + Sa^ = c„. We have at once since mi (<^„) = Sv<x etc. 

-Sv« = -Sje'r, -SViS = - cos M Sat + Syr, Svy = - sin uSat - «S/8t, 

and 

2e = FfFajS' + Vra = V^p,^ + Vra . 

Again since the double spin vector of 4>a is FVa , 

FV« = V^' f , FViS = - Far cos u + Vyr , 

FV7 = — Vat sin u — FjSr , 
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Sa^a = S^. r. S/3Vi8 = Sy^ cos u + Sar , 

SyVy = — S/3f sin u + Bar , 

2p = 2Sai. f + 2SaT, or p = Sm. T + -Sar. 

Differentiating the normal along a , we find from the Frenet formulas 

<a = Sra - SaVu, co„ = [aSa{T - V") - Sa^]fi' . 

If |8 is the normal »»„ , and y the binormal /u„ , /Sar is the torsion of the a curves. 
When dp is in the plane S^w = , da = 0, and hence for all infinitesimal 
displacements of the vertex of the trihedral of a , /3 , 7 in a plane perpendicular 
to f , a remains constant, d^ becomes ySrdp or Va^Srdp and dy is 

— pSrdp = VaySrdp . 

That is to say, the trihedral is rotated about a by the angular rotation Srdp . 
Hence for congruences orthogonal to the congruence Vdp^ = , the trihedral 
merely rotates about a by an amount equal to the projection of r on the 
tangent of the curve. In case then the f congruence is a normal congruence, 
the normal surfaces are such that displacements of the vertex on any such 
surface are accompanied by rotation about a. li dp is perpendicular to r 
there is no rotation, so that for displacements along the congruence VdpVr^ = 
there is no rotation of the trihedral at all. Particular cases would occur if 
Frr = or if r = . 

Since we can write da = — Vafi^S^dp, dfi = — V^{n^8^dp + aSrdp) , 
dy = — Vy ( /t. /Sf dp + aSrdp ) , it is evident that for any case whatever of 
displacement dp , there is a rotation given by the vector + p,^ S^dp + aSrdp . 
In case t = there is rotation about the binormal of the a curve and a is 
the common normal of the j3 and the y curves. In any case there is rotation 
about a line in the rectifying plane of a, the vector axis of rotation being a 
linear function d = p^S^ + aSr of the displacement. This function is such 
that its transverse has the normal of the a curves for a zero axis. It is also 
clear that if <t>^ is of nullity* two, then 4>^ and <(>y are of nullity unity or nullity 
two together, (since they depend directly upon the same vectors T , t ) unless 
it happen that /3' is = ;8 or = 7 in which case <py in the first instance, <f>^ in 
the second, is of nullity one greater than that of the other. 

In case (a) it is not difficult to find that the general forms of 4>p ,(f>y, </>„ are: 

<t>p = - yS^i + aSf 3 , 0v = - a8t2 + mi , <t>a= - /3Sr3 + yS^2 . 

These follow since ^^ 7 = — <t>'y^ etc. It is evident that when <p^ is given 

we also know 4>p and 4>y when further a single vector f 1 is given or found. 

Since in this case ^2, Ts are not zero and not parallel, <^„ is of nullity unity. 

* Nullity equals order minus rank, that is, here, nullity two and rank one mean the same. 
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We may have however fi either zero or parallel to ^2 or fg'. So that (f>p or <f>y 
or both may have nullity two but not nulhty three. The invariant axis of (^^ 
for the root is easily found, for the adjunct operator ^1 is 

VfiySV^z r2 = - ccSV^2 tz 

so that the invariant axis Ff 2 Ts = \^a « = 2 ^Vi V2 Saai 0:2 • This may 
also be easily found by noticing that — ^2 = ^1 7 , or ^"2 = Vi Sai 7 , and 
Ts = - Vi Sai p . 

We may have different possibilities for the other axes and roots of <f>^. 
We may have another zero root, in which case the invariant m2 = 0, that is 
SFVi V2 Vai a2 = 0, or in terms of ^2, fs, Sa^2 Ts = and a, ^2, and ^3 
are in one plane. Hence Va^2 , Va^s , V^2 f 3 are in the same direction j8" in 
the plane of /3 and 7 . 

If now the nullity is to remain unity, 0. must convert some other direction 
into that of the invariant direction ^" . This new direction cannot be in the 
plane of ^ and 7 , for if <t>^ converts this direction into /3" and annuls /3" , <(>l 
would annul the whole fi, 7 plane, and as <t>^ converts all vectors into this 
plane, ^a would annul all vectors and there would be three zero roots. Hence 
with just two zero roots and nullity unity there is a direction f not in the 
j8 , 7 plane which is converted into /3" by (f>^ . Now if ^1 f 2 is not parallel to 
j8" , Fj8" 0I f2 is such a vector. For <f>^ VjS" <j>'^ U is parallel to <f>^ V<i>'^ f 2 Va^t 
or (t>a ^2 S^2 (paa - 4>a aS^2 <f)a Tz or Ff2 F0. a0„ ^2 or finally Ff2 <x, hence 
is parallel to /3" . Now ^^ converts vectors into the plane of f 2 and fa hence 
they cannot be parallel to /3" , perpendicular to this plane. However <^^ f 2 
might vanish on account of ^2 being parallel to a. In such case fa is not 
parallel to a and will answer just as well to determine the direction in question. 
Hence we have | a unit vector in the direction of Fj8" </>! ^2 or Fj8" <t>'^ f 3 . 
This is clearly perpendicidar to 0l ^2 , and to /3" , hence is the line in the plane 
of a , f 2 , Ts which is perpendicular to <f>'^ f 2 . But ^1 a = , hence if /3' is 
the intersection of the plane of a , f 2 , Ts with the /3 , 7 plane, </>! will give for 
any vector in this plane a multiple of <^^ /3' so that if ^3' = jS cos w + 7 sin m , 
the direction f is in the plane and is perpendicular to — f 3 cos u + ^2 sin u . 
It may be written V^" ( — f 3 cos u + ^2 sin m) . Now there must under 
the hypotheses be one invariant axis in the |8, 7 plane not in the direction 
of ^" , and with a root not zero; let such direction be 8 . Then 

0^ = _ g8S8' + h^" SVfi" ^' 

where 8' is in plane of |8 , 7 and = Fa/S" , hence 

(t>^= {- gdSa + h&" 8^' ) F/3" , and <t>'^ = F/3" {gaSh - hfi'S^") , 

r2= - V^" {gocSby - h&' S^" y) , 

r3 = V&" {gocSb^ - h^' 8^" &) . 
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The operator 6 

8. It will not have escaped attention that the three operators ^„ , 4>e > <t>y 
are intimately connected and have elements in common, as for instance the 
vectors f i , ^2 , Ts • In fact it is evident at once that if we let 

e= -aSti-m2-yS^z, 
then 

<^. = - Vad, ct,, = - vpe, <Py= - Vyd, 

and that we can also write (since Sa<pc, = , etc.) 

We may arrive at ^ from another starting point however. If we consider 
the trihedral of the three unit vectors a, /3, 7 moving from one position to 
another, the original position having been oo, /3o, To, the position a, fi, y 
could be produced by a rotation q {)q, where g is a unit quaternion, and q is 
the conjugate unit quaternion, that is 

a = qooq, ^ = q^o'q, y = qyoq. 

The axis of q is the axis of rotation, and the angle of q is half the angle of 
rotation. It is important to notice that g is a definite function of a , /3 , 7 .* 
If now we displace the vertex by an amount dp , the differential change in 
the rotation would be given by 2Vdqq or 2dqq, since Tq is 1 and Sdqq = 0. 
This is a linear vector function of dp , say 6dp , that is, 

ddp = — 2Sdp^i • qiq. 

The displacement of the end of a would be V6{dp)a = da, so that 
^. = — Vad , and so for the others. From this it would follow at once that 

= a0„ + a times a linear scalar function = a<f>a — aS^i. 

We can easily now arrive at the form given above. We notice that for any 
■jj = aa + 6j8 + C7, Trj = 1 , a, b, c constant, we must have 6 = ri<f>^ — r]S^ ■ 
This prevents 6 from being arbitrary absolutely. 

9. If for q we substitute qa, we see at once that if a is a constant quaternion 

5 is not changed. That is we may use for a, ^, y any .other trihedral of 
mutually orthogonal unit vectors, which maintain constant angles with 
« , j8 , 7 in all positions. If however a is not constant then 

^1 = ^ + qcaq , where w = — 2S ( ) Vi • «i a • 

6 is the operator that converts a displacement dp into the resulting instan- 
taneous rotation, represented by an axis and a length on it which measures 

* See Joly, Manual of Quaternions, p. 26, Ex. 6. 
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the rotation rate; that is, for a displacement in the direction ri , the rotation 
has the direction and rate of rotation given by dr] . When this is multiplied 
by ds , the length of the displacement, it gives the instantaneous rotation. 
It might be called a rotation derivative. 

10. By using a fundamental identity of quaternions we have 

— ( ) Sa^a = Vaai SVi + Vai Vi Sa + FVi aSai . 

Adding similar forms for /3 and y we have 

( ) [ SaVa + SfiVP + SyVy] = a^a + )30s + T'^v " ^^Vi «-Sai 

- VVi ;8S/3i - FVi ySyi - Vai Vi Sa - Vfii Vi S^ - Vyi Vi Sy • 

But we have always ( ) = — aSa — ^Sfi — ySy so that operating on this 
with V > and transposing, 

FVi aSai + VVi /3S/3i + FVi 7S71 = - V'Vi ai Sa. - VVi jSi SjS - FVi Ti Sy . 

Hence if we set, as in § 5, 

223 = SaVa + S^V^ + SyVy 

we find from the above the important formula 

2pi) =26 + 2FVi ai Sa + 2FVi ^1 S/3 + 2FVi 7i Sy , 

or transposing and limiting V to the first following vector, 

e = p - VVaSa - FVjSSjS - VVySy , 

which gives 6 in terms of the curls of a , ^ , 7 . 
This may also be written 

e = p - VavSa - Fj3vS/3 - F7VS7 , 

where V acts on the following vector, and this reduces at once to 

p + Va(l>'a + V^<l>'^ + Vy<t>'y. 
Whence we have 

e' = p - <t>aVa - (fnV^ - <i>y Vy . 
We now have 

fi = pa + ^3 7 -0^/3, U = V? + <t>y 01. - <i>^y , ^3 = py + <t>aP -<i>ea. 

Since <^g = Vy(t>a — Va4>y, etc., these reduce to 

fi = pa + Fv« + 2Fa£, r2 =p;8 + FV|8 + 2F^6, 

Ts = P7 + Vvy + 2F7«. 

These determine the vectors fi, ^2, fs in terms of a, /3, 7 respectively, since e 
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is already known in terms oi a, fi , y . If j; preserves invariable angles with 
a,|8, y, 

dr] = pr] + VVv > 6' T] = prj + VVv + 2Vr]€ . 

11. The first scalar invariant of 6 , designated by mi is 

Sp - Sa'Va - S/3V/3 - S7V7 = P ■ 

From this we can form the function 

Xs = mi - e = FVi ax Sa + VVi |8i S^ + VVi 7i Sy . 

It is obvious at once that the curls are given in terms of 6 in the form 

Vv<x = ^ Xfl « , etc. 

The vector called the spin-vector of ^ , is the same as e in § 3, 

e« = M VaVVa + VfiVvfi + F7FV7] 

I (<^. « + <^3 ^ + </'y 7) = - H ccSVa + I3SVI3 + 7SV7] . 

A necessary and sufficient condition that e = is that the three convergences 
vanish. If the curls vanish, e = 0, and the divergences vanish, which we 
would also know from the equations /Sv« = S7V/3 — S^Vy , etc., so that the 
vanishing of the curls is a sufficient condition, for the vanishing of e . How- 
ever as p also vanishes, 6 = in this case whereas the vanishing of e means 
that 6 = 6', that is, 6 is self-transverse. Hence the vanishing of the curls 
is not a necessary condition that e = . 

When ^ = each of the operators 0„ , 4>^ , <t>^ vanishes, and a , ;8 , 7 are 
constant. When 6 is self-transverse the three sets of curves which have as 
principal normals the three vectors a, /3, 7 respectively, are such that per- 
pendicular curves of the same set have equal curvatures and one normal is 
opposite to the other; that is, taking the normal curves of a for instance, 
if one is concave in the direction of a the perpendicular curve is concave in 
the direction of — a . 

If we subtract the last form of e from the first we have 

VaVa + VfiVP + 1^777 = . 
Adding 2p we have 

q:V« + |3V/3 + 7V7 = 2p . 

12. Since 6a = VVa + pa, it is clear that if Fva = then a is an invariant 
axis of 6 , and since 0„ a = , the a curves are straight. Hence Vva = is a 
sufficient condition that the a curves be straight. But conversely if a is an 
invariant axis, either FVa = or else VaV^^a = = <t>aa hence the curva- 
ture of the a curves is zero, and they are straight lines. The necessary and 
sufficient condition is therefore F • aV'':ja = 0. If 17 is any vector which 



402 JAMES BYRNIE SHAW [October 

maintains constant angles with a , j3 , 7 then drj = FV; + pv > and if rj is an 
invariant axis for 6 at all points, so that Vridr] = = VriV^ri , then the r} con- 
gruence consists of straight lines, and displacements along them would be 
accompanied by rotation about them. Since da is the instantaneous rotation 
that accompanies a displacement along a, the component of this rotation 
along a is — aSada . In case j8 is the principal normal of the a curve and 7 
its binormal, for all positions, then — Sada is the rate of rotation of the 
osculating plane and is therefore the torsion of the a curve. However it is 
to be remembered that in this expression of the torsion 6 is dependent upon 
the normal and binormal. 6 does not change however if for the normal 
and the binormal we substitute ;8' , 7' two perpendicular unit vectors in the 
plane of j8 , 7 which maintain constant angles with them. Since — Sada = p 
— Sa\7a we see that when a is the principal normal of the /3 curves and of 
the 7 curves we have 

— Sada = p — tp — ty. 

This relation will hold if j3 and 7 are the tangents of any two perpendicular 
normal curves of a . We may write this in the form 

t$ + ty = P + Sada . 

In case a is the direction of the principal normal of a curve whose tangent 
is 18 we must have as the necessary and sufficient condition 

Va<t>^ fi = VaVm = ^Sad^ = Sad^ = S^d' a. 

Likewise if it has the direction of the principal normal of a 7 curve we must 
have Sady = 0, hence 187^' a = also. Under these conditions 

Vd' aF/37 = or Vad' a = . 

Hence the necessary and suflScient condition that a is in the direction of the 
principal normal of the j3 curves and the 7 curves, and hence of all curves 
whose tangents maintain constant angles with /3 and 7 , is 

Vad' a = . 

13. Since 2d = — aui SVi — j8j8i SVi — T7i SVi we find at once that 

2eV\7 Va2 ai SVi V2 - VV^2 /3i SVi V2 - F72 71 -S Vi V2 , 

where V on the left operates only on d and where each term on the right con- 
tains a vacant place for the operand of the linear vector operator. But this 
is the same as 

But since <t>a = — Vad we have for all X , /( , 
^; V\fi = V4>^\<S>^n = VVadWadu aSaVeXd/x a-So^i FX/t 
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SO that fa. = — aSoc^/'e with similar forms for the others. Substituting these 
we have 

^: = - ew. 

This is an important form. When it vanishes identically we have the 
necessary and sufficient condition fulfilled that 6dp be an exact differential. 
In this case we could write 6dp = — SrfpV • "^ where <r is properly chosen. 
That is, the form of d can be reduced to B = — S()V-o-. Hence ■p = 8^a , 
2e = VV<T, SVa = SaW, SV/S = SPV<r , SVy = S yVff , e„ f . = - Va4>„ a, 
etc. This particular case leads to some very interesting possibilities. 

From the equation above we have 

1^» = VVi e'l . 
From this it is easy to arrive at the second scalar invariant of 6 , for it is the 
first scalar invariant of ^«, and it is easy to prove that for the general case 
of an operator like this the first scalar invariant is — 2»SVe . 

Hence rrii = — 2Sv« • This can be verified by the direct calculation from 
the first or other form of . 

Since ^9 5 = ms ( ) , we have taking the first scalar invariant of this 

3m3 = 25 Vi isi's , 
where the subscript is removed after the differentiations. When ms = 0,6 has 
at least one zero root. In case ^e is not also zero nor m2 zero, there is but one 
zero root. The invariant axis of this root is ^9 5 where 5 is any unit vector. 
If also m2 = but \[/e is not zero, zero is a repeated root with one invariant 
axis, but S^ has a plane of invariant axes. The other root is then mi = p. 
It ms = , \l/e = , since there must be a repeated root, we must also have 
TTh = 0, and the zero root has two invariant axes, and hence a whole plane 
of invariant axes. In case all the scalar invariants are zero, there is a triple 
zero root with one invariant axis, ^ has an invariant plane, and 6^ annuls all 
vectors. If 5* = 0, there is a triple zero root, with an infinity of axes in a 
plane. li 6 = any line is an axis. 

14. Other expressions for the scalar invariants ma and ms are as follows: 

TO2 = - p2 - SaVv^VVi - S^VVyVVoc - SyVVoiVvfi, 

ms = pm2 — SFVaFVjSFVT • 
We also have 

26 + 2p = aVVa + ^Vv^ + 7FV7 , ie = VVVaa + VVvPP + VVVy^ , 
2p - 2e = VVaa + FVi^jS + FVtT . = VaV<x + VfiV^ + F7V7 , 

, , . VVaSVa. - FV/SSV/S - FV7SV7 

^(e) = pe 2 ■ = -«*« = «*«' 

= pe - |FaFFV/3FV7 - IV^VVVyVVa - ^VyVVVaVv^, 
Sede = - Seej, . 
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In all the above formulas V operates only on the first following symbol. 
Another class of formulas which are useful are the following: 

VSVa = 2vSae = - 2^1 € + 2Vi Sei a 

= 26' Vea + 2Vi Sei a = d' da - d''' a - 24>\ a, 
SaVSVa + SPvSv^ + SyVSVl = 4:6^ - 28ve, 
VavSva + F/3vSV;8 + VyvSvi 

= - 2Vae' Vae - etc. + 2Va ( Vi Sei ) a + etc. 

= - 2FVj87e' Vae - etc. + etc. 

= 266 - 2pe - 4Fv« . 

Numerous other relations can be written down easily. 

15. The linear vector function 6 has at least one invariant axis, say |i. 
This axis is such that ^^i = ^ri ^i . There is therefore at every point one 
direction at least such that displacements in that direction are accompanied 
by rotation about the same direction, though the rate gi may be zero. Con- 
sequently these directions determine at least one congruence of curves such 
that if the vertex of the trihedral travels on any one of the curves the trihedral 
will rotate about the tangent of the curve, the amount of rotation depending 
upon the position. The motion of the trihedral is a sort of screw motion. In 
case the roots of 6 are all distinct, they have shear regions of one dimension, 
and there will be at each point three curves, the whole constituting three 
congruences, such that displacement along any one of the curves is accom- 
panied by a screw motion along the curve. Since any point can be reached 
by displacements on the three congruences, any displacement and consequent 
rotations of the trihedral can be analyzed into three successive screw motions. 
If one of the distinct roots is zero, the rotation about the axis is of zero mag- 
nitude and the displacement is accompanied by a mere gliding, the trihedral 
remaining parallel to itself. If there are two zero roots which have distinct 
invariant axes, that is if they have shear regions of one dimension each, they 
may be any two distinct directions in the plane of the two shear regions. In 
such case the plane containing all the invariant axes would envelope a surface 
such that displacements from any point of the surface to any other point 
would be accompanied by no rotation whatever. If we have the similar 
case for all three roots, 6 = and any motion is a mere translation. 

When one of the roots is repeated, even a zero root, but the shear region is 
of two dimensions, so that there is only one invariant axis, then 6 converts 
any other definite direction in the shear region into a multiple of itself, which 
may be zero, plus a multiple of the unit vector of the invariant axis. Hence 
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at each point there would be one tangent line of a curve belonging to a con- 
gruence along which the motion is a screw motion on the curve, and there is 
also a tangent to a curve for which displacements are screw motions along the 
curve and also a rotation about a straight line which would be in the direction 
of the tangent to the curve at the point belonging to the first mentioned 
congruence. In case this double root is zero, there will be a congruence for 
which the trihedral simply moves parallel to itself, and there will also be a 
congruence, for which the trihedral moves parallel to itself with a superimposed 
rotation at each point about the tangent of the curve at the point which belongs 
to the first congruence. 

When there is a triple root we may have the preceding case with a third 
congruence for which the motion is a screw motion on the curve, this con- 
gruence and the invariant congruence of the same root determining a surface 
for which displacement along any curve lying on it is accompanied by a screw 
motion. The third congruence cuts through these surfaces and displace- 
ments along it would be screw motions accompanied by rotation about a single 
tangent to a curve lying on the surface. If the root is zero the rate of rotation 
is zero and the motions become translations, save for the last case. 

Finally if there is a triple root with a shear region of three dimensions, 
there is a congruence for which the motion is a screw motion, another con- 
gruence for which the motion is a screw motion accompanied by a rotation 
about a tangent to a curve of the first congruence, and a third congruence 
for which motion is a screw motion accompanied by a rotation about a tangent 
to a curve of the second congruence. In case the root is zero, the screw 
motion degenerates into a translation, but the rotations remain. 

In the cases of shear region of dimensions two or three the congruences 
accompanied by rotations are not unique, for it is clear that if 

0^2 = 01^2 + Hi , and e^i = giii, 
then 

6(^2 + xii) = gii^i + x^i) + Hi, 

and similar equations hold for a shear region of order three. 

16. If 6 has no zero roots, then ^pe a must be an axis for 0„; and similarly for 
j8 , and y , for </>„ ^« a = — Vad\p0 a = — Vaa = . In case 6 has a zero 
root the corresponding invariant axis or axes are also zero axes for each 0„, 
<l>p, 4>y. In this case if yp^ a does not vanish it is also an invariant axis for a 
zero root of all of the operators d and <!> . In case there are two zero roots 
of d , then ypo may also vanish identically, and there is a whole plane of invari- 
ant axes of 6 and each <i> , or ^^ may not vanish identically, in which case 
again ^« a is an invariant axis for <t>^ etc. The case of three zero roots for 6 
gives \(/g = 6^ , and if this does not vanish identically it gives 4'e <x an axis of 
<^„ as before. If also however 6^ = 0, and 6 7^ 0, then 6a is an axis for 0„ . 

Trans. Am. Math. Soo. 27 
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If a is an invariant axis of 5 it is a zero axis of 0„ , etc. But in this case 
Vada = , hence the curvature of curves of the a congruence is everywhere 
zero and they are straight lines. The condition Vada = is evidently the 
necessary and sufficient condition that the a congruence consist of straight 
lines. Also Vad' a = is the condition that a is everywhere the principal 
normal of the ^ and y curves. The form of 6 in this case is special, being 

e = - zaSa - ^S{u^ + vy) - yS {u' ^ + v' y) . 

Therefore Vea = and Vede = are necessary and sufficient conditions that a 
congruence of straight lines be principal normals to the orthogonal con- 
gruences ^ and y . 

17. In case Vad' a = 0, that is when a is an axis of 6' , a is the principal 
normal of the /3 and the y curves so that 

tp = 8y<f>^p = - Spep, t,= - 8^<f>^ y Sydy . 

From these we have 

tp-ty = 8y<t>^^ + S/3<^,7 = 5aFi3(<^. + <^:);8. 

Since we already have seen that 

tfi + ta = Sa^a , 
we have at once 

2tp = SaVa + SaV^{<i>a + <^l)/5, 2<, = SaV« - 8aV^{4,^ + 0:);8. 

Let us hold a constant now and vary ^ . This torsion is an extremal if 

= 8aVd^{<f>^ + <t>:)^ + SaV^{<t>^ + <t>:)dp. 

But dp is parallel to y , since a is fixed; hence 

0= -Spi<f>^ + <t>:)P + Sy{<t>^ + <j,:)y, 

Hence we have for a maximum or minimum 

/S;80„;8 = Sycf)^ y or SyvP = - S^Vy = - 5 SVa 
or 

Sydp + S/307 = . 

This may be interpreted to read: the curvatures for the curves of extreme 
torsions are both equal to — 5 the convergence of a . Hence the vectors ^ of 
extreme torsion are in the directions of equal curvature for the perpendicular 
vectors ^ , y . 

18. The curvature of a ;8 curve is easy to find when a is the normal of all 
the /3 curves, for it is 

-8a(l>pp = 8p(l>,p. 
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This is an extremal when Sdp<j)^ P + S/3</>„ dfi = 0, or since dfi is parallel to y , 

Sm = Sydy . 

That is, the direction of curves with extreme curvature is such that the torsion 
equals that of the perpendicular curve and is hence = 5 Sa^a . 
19. The value of the torsion of the ^ curve is 

2tp = SaVa + (-8/3^/3 - Sydy) = SaVa + SaV^i<t>^ + <t>'a)P 

and if we let ;8 be a direction of extreme curvature, 8^6^ — Sydy = . Sup- 
pose /3' s any other direction, and 

^3' = jS cos M + 7 sin M , 7' = a^' = — /3 sin w + 7 cos u . 
Then 

S/3' e^' - Sy' dy' = 2S (18^7 + 7^18) sin 2m. 

Hence for any torsion 

2tp, = SaVa + 25 (18^7 + yd^) sin 2m. 

Again the value of the curvature of the /3 curve is Syd^, and if we measure 
u now from the direction of extreme torsion, we have 



or 



-S7' e^' = Syd^ - HS^ep - SyOy) sin 2m 
2S7' e^' = -Sv« - {Sm - Sydy) sin 2m 



as the value of the curvature of the 0' curve. 

It is easy to see if we set /3' = |8 + 7 , 7' = /3 — 7 , that when 

Sm - -87^7 = 
we have 

Sp' dy' + Sy' d^' = 0, 
and when 

Spey + Syd^ = , S/3' dfi' - Sy' ^7' = . 

Hence the two sets of extremal directions bisect each other. 

20. Recurring to the form 6 = q:0„ — aS^ we need to notice that this form 
alone with f purely arbitrary will not give us d , since it is too general. The 
derivation of 6 shows that a similar form must hold for any unit vector -q , 
instead of a , if 77 is fixed with reference to the trihedral of a , /3 , 7 . Hence J" 
is not purely arbitrary. In fact, we have from the form given in § 12 for i^# , 
when 6 is taken in the form just above, 

eVv = - Fa2 «i SVVx V2 + ai SVi r ( ) + «SFvf ( ) = - t^i, 

whence 

^<, = 2^„- V^<t>:-Vv^Sa. 
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But forming xf/^ FX/x = Vd' \d' n directly we get 

It follows that we ought to have Fvf = —</'««• This we can prove directly 
as follows : 

V\7^ = - Fv</)^ T = FVV2 Sp2 7 = VVi V2 S^^ 71 = VVi <t>'p 7i • 

Now we have from § 10, (^3 = — <f>'^Vy + (f>'^Va, and substituting we have: 

Fvr = VVi [ <f>: F71 Va^ - </,; F71 F;87] 

= VVi [ <t>'a ^Sayi - 0; aSpyi - 4>', yS^ji + <j>'^ pSyyi ] 

= — FVi <t>a PSai 7 = V<t>'a 7<^1 /3 = — ^„ a . 

The form above gives simple forms for the invariants of 6 . We have 

mi = p = SaVa - /Saf , ma = SaVt + S^Va, rris ^ S^Vf , 

2e = a/SVa + <^. a + Faf , 2^6 = (A, VVa - ^1 f - fSVa . 

The discussion above can now be stated in terms of a and f . For instance : 
6 cannot have a zero root unless rtiz = , that is SfVf = 0, and the f con- 
gruence is normal to a set of surfaces. We need not elaborate these results 
further however. 



